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I.  INTRODUCTION 


A.  OBJECTIVES 

Finite  state  space  semi-Markov  models  find  application  in  a  variety  of  areas  such 
as  queueing  theory,  reliability,  and  clinicai  trials.  The  application  of  these  models  often 
centers  on  the  distribution  of  a  first-passage  time  to  a  state  or  a  set  of  states 
representing  for  example  the  lifetime  of  a  system  or  marking  the  end  of  a  busy  period 
for  a  server.  Suppjose  that  the  observations  of  the  path  of  the  semi-Markov  process  are 
all  that  is  known  about  the  process. 

This  thesis  reports  the  results  of  a  simulation  experiment  to  compare  various 
parametric  and  nonparametric  estimators  of  the  natural  logarithm  of  the  probability  a 
semi-Markov  process  does  not  enter  a  particular  state  before  a  given  fixed  time  t.  In 
what  follows  we  wail  use  "In  probability"  for  "  natural  logarithm  of  the  probability 
The  specific  semi- Markov  model  and  estimators  considered  are  given  in  Chapter  II. 
Chapter  III  contains  the  details  of  the  simulation  experiment  and  results.  Conclusions 
from  the  study  are  given  in  Chapter  IV. 

B.  SCOPE  OF  THE  THESIS 

The  purpose  of  this  thesis  is  to  use  simulation  to  compare  estimators  and  their 
confidence  intervals  for  the  In  probability  a  particular  semi-Markov  process  does  not 
enter  a  particular  state  before  a  given  fixed  time  t.  The  particular  semi-Markov  model 
and  estimators  considered  are  given  in  Chapter  II.  A  simulation  study  comparing  bias 
and  standard  errors  for  these  estimators  was  reported  in  Gallagher  (198b)  [Ref.  1].  In 
this  thesis,  we  are  primarily  interested  in  comparing  confidence  interval  procedures. 

An  estimation  procedure  which  uses  the  least  information  about  the  semi- 
Markov  process  uses  the  empirical  distribution  of  the  observed  first  passage  times.  An 
estimation  procedure  which  uses  the  most  information  is  to  assume  a  parametric  form 
for  the  sojourn  time  distribution  in  each  state  and  a  transition  matrix  to  describe  the 
transition  between  states:  the  parameters  of  the  sojourn  time  distribution  and 
transition  probabilities  can  be  estimated  using  maximum  likelihood.  An  estimation 
procedure  requiring  less  information  uses  nonparametric  estimators  of  the  sojourn  time 
distributions  and  the  maximum  likelihood  estimators  for  the  'ransition  probabilities. 


10 


In  many  cases,  parametric  assumptions  concerning  the  sojourn  time  distribution 
are  difficult  to  justify.  It  was  demonstrated  in  Park  (1986)  (Ref.  2]  that  incorrect 
parametric  assumptions  may  lead  to  very  biased  estimators.  Hence,  a  nonparametric 
estimation  procedure  may  be  preferred  to  a  parametric  one  when  actual  data  is  used. 
However,  the  nonparametric  procedure  can  be  expected  to  be  less  efficient  than  a 
parametric  one  provided  the  parametric  model  is  correct. 

The  thesis  is  organized  as  follows.  In  Chapter  II,  the  nature  of  the  problem  is 
described  and  several  parametric  and  nonparametric  confidence  interval  estimation 
procedure  are  introduced.  In  Chapter  III,  the  simulation  experiment  is  described  and 
results  are  given.  Conclusions  drawn  from  the  simulation  study  are  given  in  Chapter 
IV. 


II.  NATURE  OF  THE  PROBLEM 


A.  DESCRIPTION  OF  PROBLEM 

The  semi-Markov  process  model  used  in  the  simulations  to  compare  the 
estimators  is  as  follows. 

Suppose  we  observe  N  individuals.  Let  Xt  (i)  be  the  state  of  the  i-th  individual 

at  time  t.  We  will  assume  (Xt  (i),  t£  0}  i  =  l,  2 .  N,  are  independent  semi- 

Markov  processes  with  three  states  (  0,  1,  2  }  having  the  same  distribution  as  {Xt  .  t^ 
0}.  All  individuals  start  at  t  =  0  in  state  1,  (XQ  =  1).  An  individual  stays  in  state  1  for 
a  random  length  of  time  having  distribution  function  F,.  Upon  leaving  state  1.  the 
process  transitions  to  state  0  with  probability  0  and  to  state  2  with  probability  1-0  . 
If  the  process  transitions  to  state  2,  it  spends  a  random  length  of  time  there  having 
distribution  function  F,.  From  state  2,  the  process  transitions  to  state  1  with 
probability  1.  State  0  is  an  absorbing  state.  Once  in  it.  the  individual  never  lea\es. 

For  all  individuals,  the  entire  path  of  transitions  and  sojourn  times  are  observed 
until  the  time  of  absorption  in  state  0  (Fig.  2.1). 

Let 


D  =  inf  (t£  0,  Xt  =0}, 

be  the  entrance  time  to  state  0  (or  the  time  of  death). 

The  problem  is  to  estimate  the  logarithm  of  the  first  passage  time  survivor 
function  P  (D>t)  for  fixed  time  t  from  the  data  obtained  by  observing  the  N 
individuals. 


B.  THEORETICAL  DEVELOPMENTS 

Two  estimators  for  P(D>t}  will  be  described  in  this  section.  They  are  the 
maximum  likelihood  estimator  and  the  asymptotic  renewal  estimator  from  a  paper  by 
P.A.  Jacobs.  [Ref.  1:  p.2-6] 
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Figure  2.1  Transition  State  Step  (Sojourn  Time). 

1.  Maximum  Likelihood  Estimate  for  Continuous  Time  Markov  Chain 

In  this  subsection,  the  maximum  likelihood  estimator  for  P|D>  ti  will  be  given 
for  the  special  case  in  which  the  sojourn  time  in  state  i  is  exponentially  distributed 
with  mean  l.p^  (i*  I.  2). 

Let  R  be  the  number  of  transitions  from  1  to  2  for  one  individual.  The  log 
likelihood  function  for  the  individual  is 

L  *  R  In  (1-0)  +  In  0  +  R  In  p,  (eqn  2.1) 

-i-  (  I  +  R  )  In  Pj  -  PjTj  -  p,T2 

where  Tj  (i*  1,2)  is  the  total  time  spent  in  state  i  before  death  . 

The  maximum  likelihood  estimators  are 

$  a  - ! -  (eqn  2.2) 

1  +  R 

A  1  ■+■  R 

P[  *  -  (eqn  2.3) 

T1 

13 


V«7V. 


(eqn  2.4) 


A  R 
Pj’“ 

Further  R  is  geometric  {8}  so  E{R}  =  (1  -  0)/0.  Since  the  N  individuals  are 
independent,  the  maximum  likelihood  estimators  using  the  data  for  all  N  individuals 
are 


N  +  R 


(eqn  2.5) 


a  N  +  R 

Pi  _  — - -  (eqn  2  6) 

T1 


A  R 

P2  *  — r-  (eqn  2.7) 

T, 

* 


r-*  .-*✓ 

where  R  is  the  total  number  of  transitions  from  1  to  2  for  all  N  individuals  and  T;  is 
the  total  time  spent  in  state  i  for  the  N  individuals.  For  the  estimators  based  on  data 
for  all  N  individuals. 


Var  (0}  = 

I  t0)_I 

(eqn  2.8) 

A 

Var  (Pj)  = 

I  (P!)"1 

(eqn  2.9) 

A 

var  fp2}  - 

I  (P:)”1 

(eqn  2.10) 

where 


1(0)  =  N  E  f 


R 

( 1-0)* 


(eqn  2. 1 1) 
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KPP  -  N  E  {■ 


1  +  R 


(eqn  2.12) 


Pi* 


■} 


I(P2) 


N  E  (- 


P2 


(eqn  2.13) 


with  N  being  the  number  of  individuals. 

Let  D  be  the  time  of  entrance  to  state  0  for  an  individual.  Fix  t  and  put 


S  -  S(t)  =  P  {D>  t}; 


then 


S(t) 


X,+  p,  ,  .  X,  +  p,  Op, 

=  {  — 2-.-^  EXP  (X2  t)  -  EXP(Xj  t)  }  p- 

A,,  A. 


xrx2 


(eqn  2.14) 


where  X,,  X-,  are  the  roots  of  the  equation 


GpiP:  +  y  (Pi  +  P2)  +  T  =  0  • 


(eqn  2.15) 


A  parametric  estimator  for  S(t),  the  survival  function,  is 


P-, 


A  A  A  A  AA 

A  A-,+  P-,  ,  X,  +  p,  0p« 

PM  [E»  t}  =  {  a--  EXP(X2t)  -  EXP(Xjt)  j- 

Xj  Xj-X2 


A 

X, 


(eqn  2.16) 


where  Xj  and  X2  are  roots  of  the  equation 


AAA 


0PiP:  +  y  (pj  +  p2)  +  y ■  o  • 


(eqn  2.17) 


Since  the  maximum  likelihood  estimators  are  uncorrelated,  the  asymptotic 
a  a 

variance  of  P-yj  (D>  t],  which  will  be  denoted  by  Var^j.  is  approximately 


VarM  !PV1  'D>t)|0.p1.p2; 


(eqn  2.18) 
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3S  3A1 


3S  3X2 


3a  30 


3 a?  30 


(eqn  2.19) 


3 S  _  S_  33_  aX1  +  3S_  ^2 

3o1  3i  3Ai  3p)  *  3A2  3p1 


(eqn  2.20) 


i£_  .  _ - ] — -  _L_eAlw  .*££_  .  (eqn  2.21 

lAi  A2-  a2  A  j  *  3 a ^  op^  3 A2  3 *5 


V  A  ,  A  } 
i  2 


e°1  r»2  A,-. 

v  A,  - a„ j  L ,  2 
J  C  A, 


-J  aA'-  -  -  1  )te  ‘"  ]  (eqn  2.22) 


V*2 


r^TT  [*-4  «*'*  *  (-T  *  ’)«**'•]  "qn::•:3, 

A  1  A  '  y  C*  A  ^ 


1  2  X 


-0.  P, 

I  ' 


(eqn  2.24) 


30  p1*p2*2a1 


3^i  <3p2 


(eqn  2.25) 


«P,  P,*P2*2>i 


3^i  _*i* 

7TZ  *  P,-P2*2i! 


(eqn  2.26) 


^  A  A 

The  asymptotic  variance  of  In  Pyj  (D>  t}  »  Varyj  /  (P\j  {D  >  t> )*. 

2.  Asymptotic  Renewal  Estimator 

This  subsection  describes  an  estimation  procedure  for  P  {D>  t},  which  is  based 
on  an  exponential  approximation  to  the  survivor  function  P{D>t}.  Details  of  the 
approximation  which  is  obtained  by  asymptotic  renewal  theoretic  results  is  given  in 
Jacobs  [Ref.  3].  The  approximation  improves  as  t  becomes  larger. 

A 

Let  0  =  N  (N  +  R  )  be  the  maximum  likelihood  estimator  of  0  (equation 

A 

2.5).  Let  ipj(a)  be  the  empirical  transform  of  the  sojourn  time  in  state  i  ;  that  is.  if 
S;(i),  S,(i) . Syj  (i)  denote  ail  the  sojourn  times  in  state  i  for  the  N  individuals,  then 

t 

A  M‘ 

<p-  (a)  =  (1/M;)  V  EXP  (a  x  Sj.  (i)  ).  (eqn  2.2'") 

k  =  1 


The  asymptotic  renewal  estimator  of  the  survival  distribution  |Ref.  3]  is 


A  AAA 

Pa(D>  t}  =  (b  p)  EXP(-K  x  t). 


where  K  is  the  solution  to  the  equation 


A  AAA 

Ra)  *  (1-0)  q>j(a)  <p;  (a)  -  1  =  0  . 


n  M.  M, 

A  (1-0)  _l  -  A 

"  J,-  -  ;Sk<»+S|t2*»  EXPfU^n+SjCW. 

1  ‘  :k=  ij=  i 


A  A 
(pj  (K) 


M, 


(1/M;)  V  EXP  (K  x  Sk(l)  ). 
k=  i 


and 


A 

b 


-A 


K)  <Pj(K). 
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(eqn  2.28) 


(eqn  2.29) 


(eqn  2.30) 


(eqn  2.31 ) 


(eqn  2.32) 


Note  that  when  a*0,  the  left  hand  side  of  the  equation  (2.29  >  is  -  8<  0.  As 
a  increases,  the  left  hand  side  of  the  equation  increases  to  sc.  thus  there  is  a  unique 

A 

solution  K.  The  solution  may  be  found  by  numerical  search.  One  possible  numerical 
search  procedure  is  the  golden  section  search  method. 

3.  Empirical  Distribution  estimate 

Let  dp  d, . d^-  denote  the  observed  times  of  absorption  in  state  0  for  the  \ 

individuals.  A  binomial  estimator  for  the  survivor  function  is 

\ 

PBiD>t}  -  (1/N)V  i(tt0o)(d.)  ieqn  2.33) 

k=  t 

where  l(t,x)  (<*•,)*  1  if  d-x >  t, 

0  otherwise. 

C.  CONFIDENCE  INTERVAL  PROCEDURE 

Suppose  X  is  a  random  variable  whose  probability  law  depends  on  an  unknown 

parameter  0.  Given  a  random  sample  of  X:  Xj,  x-> .  xn  the  two  statistics  lower  ( L > 

and  upper  (LT  form  a  L00(  l  -  ct  )%  confidence  interval  for  0  if 

PlLSfl^  U)  £  1  •  a. 

Procedures  to  obtain  confidence  intervals  for  the  point  estimates  for  In  P(D>  tl 
are  given  below. 

1.  Maximum  Likelihood  Estimator 

The  maximum  likelihood  estimator  is  asymptotically  normal  as  the  number  of 
individuals  X  becomes  large,  [Ref.  4].  For  a  fixed  t.  the  maximum  likelihood 
confidence  interval  for  In  Pl  D>  t|  is 

{  L  .  U  }  =  In  PVj(D>tJ  ±Zpa  ,  (v  Var^j )  (P^  !D<t;).  (eqn  2.34) 

A 

where  Var^j  is  obtained  by  using  the  maximum  likelihood  estimator  of  0.  Pp  and  p,  in 
equation  12. 11)  -  12.26).  Confidence  limits  that  are  larger  than  0  are  set  equal  to  <>. 


IS 


2.  Empirical  Distribution  of  First  Passage  Time  (  Binomial  and  Normal  Approx. 
Confidence  Interval  Methods  ) 

a.  Binomial  Confidence  Interval  ^ 

Since  the  N  individuals  are  assumed  independent,  the  estimator  Pg{D>  t} 

of  equation  (2.33)  has  a  binomial  distribution  with  N  trials  and  probability  of  success 

P(D>  t}.  The  IMSL  routine  BELBIN  was  used  to  obtain  binomial  confidence  intervals 

for  P(D>t}.  For  a  description  of  the  procedure  see  [Ref.  5:  p.390  -  391).  The  binomial 

confidence  interval  for  In  P(D>t)  is  obtained  bv  taking  logarithms  of  the  upper  and 

* 

lower  confidence  limits.  If  the  lower  confidence  limits  for  Pg  {D>  t}  is  equal  to  0,  it  is 
set  equal  to  0.000 1 . 

b.  Approximate  Normal  Confidence  Interval. 

If  the  number  of  individuals  N  is  large,  the  binomial  confidence  interval 
procedure  can  be  approximated  by  a  normal  confidence  interval  procedure  as  follows, 
[Ref.  6:  p.99  -  100  .  p.954  -  955].  The  interval  is 


{  L  .  U  }  =  Pg[D>t)  ±z1,a/:v'Varg  (eqn  2.35) 

where 

Varg  -  (  I  ••  N)(  Pg(D>t}  x(  l  -  Pg(D>t}  )} 

This  interval  {  L  ,  U  )  is  an  approximate  100(1  -  a)%  confidence  interval  for  P[D>t). 
The  interval  for  In  P(D>  t}  is  {  In  L  ,  in  U  }.  If  either  L  or  U  is  less  than  0,  it  is  set 
equal  to  0.0001.  If  either  L  or  U  is  greater  than  I.  it  is  set  equal  to  1. 

3.  Asymptotic  Renewal  Procedure 

In  this  subsection,  two  nonparametric  methods  are  described  for  obtaining 

A 

confidence  interval  using  the  estimator  P^  (D>t}.  One  is  the  jackknife.  The  other  is 
the  bootstrap.  They  are  described  below, 
a.  Jackknife  Estimation  Method 

The  jackknife  was  first  introduced  by  Quenouille  (1956)  [Ref.  7],  for  the 
purpose  of  reducing  the  estimate  bias,  and  the  procedure  was  later  utilized  by  Tukey 
(1958)  [Ref.  8],  to  develop  a  general  method  for  obtaining  approximate  confidence 
intervals. 

The  basic  idea  of  the  jackknife  estimation  method  is  to  assess  the  effect  of  each 
of  the  groups  into  which  the  data  have  been  divided,  not  by  the  results  for  that 
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group  alone,  but  rather  through  the  effect  upon  the  bodv  of  data  that  results 
from  omitting  that  srouo.  The  two  bases  of  the  jackknife  are  that  we  make  the 
desired  calculation  for  all  the  data,  and  then,  after  dividing  the  data  into  groups, 
we  make  the  calculations  for  each  of  the  slightlv  reduced  Bodies  of  data  obtained 
bv  leaving  out  just  one  of  the  groups.  [Ref.Y] 


The  jackknife  procedure  is  as  follows: 

A 

(1)  Let  Yajj  (t)  be  the  estimate  P^{D>  t}  computed  using  all  the  data. 

(2)  Let  Yj  (t)  be  the  computed  statistic  using  that  data  which  omits  the 

subgroup  where  j  =  1,  2,  ,  k.  In  this  thesis,  for  the  case  in  which  the 

number  of  individuals  N  equals  20  the  j^  subgroup  consists  of  all  data 
corresponding  to  the  j1*1  individual.  For  the  case  in  which  the  number  of 
individuals  N  equals  50  the  first  subgroup  consists  of  all  data  corresponding  to 
the  first  5  individuals,  the  second  subgroup  of  the  second  5  individuals  etc. 
Some  cases  for  N  =  50  were  run  with  the  j^  subgroup  consisting  of  the  j1*1 
individual.  The  resulting  confidence  intervals  differed  little  from  those,  obtained 
by  leaving  out  5  individuals  at  time.  Computational  considerations  lead  us  to 
use  the  fewer  number  of  subgroups  in  this  case. 

f  h 

(3)  Define  the  jin  pseudo- value  by 

Y*j(t)  =  k  In  Yalj(t)  -  (k-l)  In  Yj(t).  (eqn  2.361 

(4)  The  jackknife  estimator  Y*(t)  is 

Y*(t)=  ( 1. k)  1 Y*  (t)  +Y*  (t)  + . +  Y*  (t)}.  (eqn  2.37) 

12k 

(5)  The  jackknife  estimator  of  the  variance  of  Y#  (t)  is 

k 

S*:(t)  =  (I  (k  (k- 1 ))}  V  (Y*j(t)  -  Y*(t))2.  (eqn  2.38) 

j=l 

Tukey(1958)  proposes  that  in  a  wide  variety  of  problems  the  k  estimated 
pseudo-values  can  be  treated  as  approximately  independent  and  identically 
distributed  random  variables  [Ref.  8).  to  obtain  the  following  confidence 
interval  procedure. 

(6)  The  jackknife  confidence  interval  is  computed  as  follows. 


(  L.U}  -  Y.(t)  ±  t1<2VS*'(t) 


(eqn  2.39) 


where  tpa,2  is  the  upper  1-  a/2  critical  point  of  the  t-distribution  with  k-1 
degrees  of  freedom.  The  confidence  interval  given  by  equation  (2.39)  is  a 
function  of  the  estimated  variance.  If  either  confidence  limit  is  greater  than  0. 
it  is  set  equal  to  0. 
b.  Bootstrap  Estimation  Method 

Efron(  19791  introduced  the  bootstrap  method  for  estimating  the 
distribution  of  a  statistic  computed  from  observations,  [Ref.  10).  In  this  thesis,  the 
bootstrap  was  implemented  as  follows. 

Suppose  data  are  gathered  for  N  individuals.  Let  R[2(n)  be  the  number  of 
transitions  from  state  1  to  state  2  for  individual  n.  Let  (Sj.(i)}  be  the  collection  of  ail 
sojourn  times  in  state  i  for  all  individuals.  A  bootstrap  replication  is  generated  as 
follows:  To  generate  data  for  one  individual,  one  observation  is  drawn  at  random  with 
replacement  from  (Rp(n)};  call  the  observation  rp;  rp+  l  observations  are  drawn  at 
random  with  replacement  from  the  collection  of  state  1  sojourn  times  [Sj.(l)};  rp 
observations  are  drawn  at  random  with  replacement  from  the  collection  of  sojourn 
times  in  state  2,  (Sj,(2)}.  This  procedure  is  replicated  N  times  to  generate  bootstrap 
data  for  N  individuals.  The  estimator  In  P^  (D>t[  is  computed  using  the  generated 
data.  This  completes  one  bootstrap  replication;  B  bootstrap  replications  are  done. 

A 

The  B  estimates  of  In  PA  (D>t}are  ordered.  A  100*(l-a)%  confidence  interval  is 
constructed  using  the  a/2  and  1  -  a/2  quantiles  of  the  bootstrap  estimates  of  In  PA 
(D>t|.  If  either  confidence  limit  is  larger  than  0.  it  is  set  equai  to  0.  In  the 
simulations  the  number  of  bootstrap  replications  is  B  =  100. 


III.  ANALYSIS  OF  THE  CONFIDENCE  INTERVAL  PROCEDURES 


A.  DESCRIPTION  OF  SIMULATION 

A  Fortran  program  is  written  to  simulate  the  semi-Markov  process  and  compute 
the  confidence  interval.  All  simulations  are  earned  out  on  an  I B  VI  3033AP  Computer 
at  the  Naval  Postgraduate  School  using  the  LLRANDOM  II  Random  number 
generating  package  [Ref.  13]. 

The  data  for  the  simulated  experiments  are  generated  as  follows: 

(1)  An  individual  starts  in  state  1  at  time  0. 

(2)  A  random  number  with  distribution  FL  is  generated  for  the  sojourn  time  in 
state  l. 

(3)  A  uniform  random  number  is  generated. 

If  it  less  than  0.  the  process  transitions  to  state  0  and  the  data  for  one  individual  is 
complete.  If  it  is  greater  than  0.  the  process  transitions  to  state  2.  A  random  number 
having  distribution  function  F,  is  generated.  The  procedure  then  returns  to  step  2. 
Data  are  generated  for  N  individuals  and  are  collected  as  follows: 

( 1)  The  passage  times  to  state  0  for  each  individual: 

(2)  The  sojourn  times  in  states  1  and  2  for  each  individual: 

(3)  The  number  of  transitions  from  state  i  to  state  2  for  each  individual. 

For  each  replication  of  the  simulation  the  estimates  and  confidence  interval  in 
Chapter  II  are  computed  for  In  P[D>t).  Each  simulation  is  replicated  300  times.  The 
true  In  PfD>t}  is  computed.  For  each  procedure,  the  number  of  confidence  intervals 
covering  the  true  value  is  recorded. 

The  number  of  individuals  that  are  too  low  (  true  In  P!D>t}>L';  the  upper 
endpoint  of  the  interval  )  and  too  high  (  true  In  P(D>  t!  <  L;  the  lower  endpoint  of  the 
interval  )  are  also  recorded.  The  average  length  of  the  confidence  interval  is  computed 
as  well  as  the  standard  deviation  of  the  lengths.  The  star.  Jar  i  deviation  is  computed  by 
subtracting  the  mean  length  from  each  length,  squaring  the  -esuits.  summing  over  the 
300  lengths  and  finally  divided  by  299. 


B.  SIMULATION  RESULTS  FOR  THE  EXPONENTIAL  MODEL. 

In  this  section  results  will  be  reported  for  a  simulation  experiment  in  which  the 
sojourn  time  distributions  in  both  states  are  exponential;  Fj  is  exponential  with  mean 
l.'Pj  and  F,  is  exponential  with  mean  1/  p2-  Some  true  values  for  P{D>  t}  for  this  case 
can  be  found  in  Appendix  A. 

For  each  replication  nominal  80%,  and  90  %  confidence  intervals  for  In  P{D>  t} 
for  various  values  for  t  are  computed  using  each  procedure  of  Chapter  II.  The  times 
considered  are  t  =0.5,  1.0,  1.5,  2.0.  3.0,  4.0. 

The  confidence  interval  procedures  for  In  P(D>t)  are  the  binomial  confidence 
interval  (=  BIN  ),  and  it's  approximating  normal  confidence  interval  (=  NOR  ),  for 
the  In  fraction  of  individuals  who  have  not  entered  state  0  by  time  t;  the  maximum 
likelihood  confidence  interval  (=  MLE  ),  jackknife  confidence  interval  for  the 
asymptotic  renewal  estimator  (=  J.K  ),  and  bootstrap  confidence  interval  for  the 
asymptotic  renewal  estimator  (=  B.T  •).  For  each  procedure  the  number  of  intervals, 
covering  the  true  value  of  In  P  (D>t)  is  recorded  as  well  as  the  number  of  intervals 
that  are  too  high  or  too  low.  These  results  are  reported  in  Tables  1,  3.  5  and  7.  Next  to 
each  coverage  count  are  given  the  corresponding  coverage  proportion  in  parenthesis. 
Below  the  coverage  numbers  is  given  a  95%  confidence  intervals  for  the  coverage  rates. 
This  interval  is  computed  as  follows.  If  P  is  the  proportion  of  (l-a)%  intervals  that 
cover  the  true  value  then  a  95%  confidence  interval  for  the  coverage  rate  is 

{  CL  ,  CU  }  =  P  ±  1.96V  P(l-  P  )/300. 

If  a  (l-a  )%  confidence  interval  procedure  is  performing  well,  then  this  interval  should 
cover  about  (!-</)%  of  the  time.  In  addition,  if  a  80%  (respectively  90%)  confidence 
interval  procedure  is  working  well,  then  out  of  300  replications  between  226 
(respectively  260)  and  254  (respectively  280)  confidence  intervals  should  cover  the  true 
value.  Simulations  for  which  the  numbers  of  intervals  that  cover  are  outside  of  these 
bounds  are  given  in  bold  type.  In  the  first  column  of  Tables  1,  3,  5  and  7  the  true 
PfD>  t}  is  given  in  parentheses. 

In  Tables  2,  4,  6  and  8  are  recorded  the  average  length  of  the  confidence  intervals 
for  In  P(D>t}.  The  estimated  standard  deviation  of  the  length  is  below  the  average 
length  in  parenthesis.  If  a  procedure  is  performing  well,  it  should  not  only  have  the 
correct  coverage  rate  but  also  a  small  average  length.  The  simulation  results  recorded 
in  Tables  1  -  4  are  for  a  simulation  whose  parameter  values  are  Pj  =  1.0.  p,=  10.0  and 


0-0.5.  The  number  of  individuals  N  is  set  at  20  and  50,  representing  a  low  and 
moderate  number  of  individuals. 

Coverage  results  for  N  -  50  individuals  are  presented  in  Table  1.  The  binomial 
confidence  interval  procedure  tends  to  overcover.  The  maximum  likelihood  procedure 
has  the  right  coverage  for  all  cases.  The  two  confidence  interval  procedures  using  the 
asymptotic  renewal  estimators  undercover  for  t^l.O  and  have  the  right  coverage  for 
t>  1.5.  When  a  jackknife  confidence  interval  does  not  cover  it  is  often  because  it  is  too 
high  (true  value  <  L  ). 

Sample  means  and  variance  of  the  confidence  interval  lengths,  for  N  =  50 
individuals,  are  reported  in  Table  2.  The  average  lengths  for  the  maximum  likelihood 
and  bootstrap  are  very  close  for  t  larger  than  1.0.  The  bootstrap  and  jackknife 
procedures  have  larger  average  confidence  interval  lengths  than  those  for  the  maximum 
likelihood  estimator  but  smaller  than  for  the  binomial  procedure. 

Results  for  a  simulation  with  X  -  20  individuals  are  given  in  Tables  3  and  J. 
They  are  similiar  to  those  in  Tables  1  and  2.  However,  the  average  length  of  the 
intervals  is  larger  in  Table  4  than  in  Table  2  reflecting  the  smaller  number  of 
individuals. 

Results  of  a  simulation  with  the  same  parameters  as  the  first  but  with  the  0  = 
2/3  are  given  in  Tables  5  -  8.  Results  for  N  =  50  individuals  are  given  in  Tables  5  and 
6.  The  coverage  results  are  in  Table  5.  The  binomial  procedure  has  better  coverage 
than  in  Table  1.  The  procedure  based  on  the  asymptotic  renewal  estimators  have  the 
correct  coverage  only  for  t2:  2.0.  The  average  length  of  the  intervals  in  Table  6  are  in 
general  longer  than  the  lengths  in  Table  2. 

Results  for  a  simulation  with  0=  2,  3  and  N  =  20  individuals  are  given  in  Table 
7  and  8.  The  binomial  confidence  interval  tends  to  overcover  when  the  true  probability 
is  greater  than  0.5241  or  less  than  0.2753.  Once  again  the  confidence  intervals  based 
on  the  asymptotic  renewal  procedure  undercover  for  t  <  2.0.  The  maximum  likelihood 
confidence  intervals  give  the  correct  coverage. 


TABLE  1 


time  C.I  (%) 

coverage 

BIN 

NOR 

MLE 

j.K 

B.T 

0.5  80  % 

too  hiah 

171.06) 

381.13) 

271.09) 

1191.40) 

1511.50) 

cover 

« 

2391. SO) 

2451.  S3) 

1671.56) 

1491.50) 

too  low 

231.08) 

281.09) 

141.05) 

()(  .00) 

int.val 

(.S3;.91) 

(.75. .84) 

t  .”..86) 

( .50..0D 

(.44, .55) 

( .  7S66) 

too  hiah 

51.02) 

171.06) 

141.05) 

831.23) 

1101.37) 

cover 

2881.96) 

2T6c'.92) 

2741.01) 

2091.70) 

1901.63) 

too  low 

'(.02) 

~1.02) 

121.04) 

S(.03j 

0  .00) 
i.5S..o9> 

int.val 

(.94.. 98) 

(.89. .95) 

(.88.  .95) 

( .04.. '5 1 

1.0  SO  % 

too  hiah 

281.09) 

2S(.09) 

271.09) 

'01.23) 

SOI. 2") 

cover 

2451.82) 

245).  82) 

2451.82) 

2161.72) 

2131.71) 

too  low 

271.09) 

2 '1.09) 

281.09) 

141.05) 

"(.02) 

int.val 

1.77..86) 

( .'7..S6) 

(.77..S6) 

( .67..  77) 

(.66.  .76) 

(.6203) 

90  % 

too  hiah 

Sf.03) 

121.04) 

141.05) 

451.15) 

451 .15) 

cover 

2821.94) 

2751.92) 

2'41.91) 

2471.82) 

2501.83) 

too  low 

101.03) 

131.04) 

121.04) 

Sl.03 1 

5)  .02) 

int.val 

(.91. .97) 

(.89.  .95) 

(.88.  .95) 

(.7  S..S7) 

1 .79. .88) 

1.5  SO  % 

too  hiah 

19(  .1)6 1 

23(o)8) 

2't.oO) 

531. 1~) 

46(  .15) 

cover 

2621.87) 

2431.81) 

2451.82) 

235<‘."S) 

23  31. '8) 

too  iow 

191.06) 

341.11) 

281.09) 

121.04) 

211.07) 

(.4891) 

int.val 

1.84;.91) 

(.77;.  85) 

(.77. .86) 

I.74..83) 

(.'3. .82) 

90  % 

too  hiah 

111.04) 

2781.93) 

111.04) 

141.05) 

321.11) 

2621.87) 

231.08) 

! 

1 

cover 

2781.93) 

2741.91) 

2661.89) 

too  low 

111.04) 

111.04) 

121.04) 

0(  .02) 

1  11.04) 

int.val 

!  .89. .96) 

i  .90. .93) 

(.8S..95) 

(.84..  91 ) 

(.85. .92) 

1  2.0  son  !  too  hiah  i  181.06)  31(.10>  271.09)  471.16)  311.10) 

i  !  cover  |  260(.87)  2331.78)  2451.82)  24<)(.S0)  243(.8h 

I  I  too  low  |  22c  .07)  36(.12)  2S1.09)  13(  04)  26  .09) 

I  !  int.val  (.83. .91)  (.73..S2)  (.77..S6)  1.75..S5)  (."..85) 

'(.385.)  ! 

90  °o  :  too  hiah  9(.o3)  9(.03)  14(.05)  22(.07)  151.05) 

cover  2771.92)  2691.90)  274.91)  2"0!.9«)  2681.89) 

I  !  too  low  j  !41.0.?»  22(.07)  121.04)  Sl.03)  17(.06) 

I  int.val  •  (.S9..951  (.86. .93)  I.SS..95)  (.87. .93)  (.86.  .93) 


too  man 

201.07) 

201.07) 

2  1.00) 

36(  .12) 

281.08) 

cover 

2621.87) 

2471.82) 

2451.82) 

2501.83) 

2401.80) 

too  low 

!  181.06) 

331.11) 

281.00) 

141.05) 

351.12) 

int.val 

!  (.84. .91) 

(.'8. .87) 

(.77. .86) 

(.79. .88) 

(.75. .85) 

too  hieh 

1  13(.o4) 

13(9)4) 

141.05) 

20(.()7) 

1 1(  .04) 

, .  v  vj  r 

2831.94) 

2601.90) 

2751.92) 

2741.91) 

26'(.S9) 

roo  iow 

41.01) 

1 81.  On) 

1  11.04) 

0<  9)2) 

2219)7) 

:nc.  ...1 

(.9  2. .97) 

(.86..  93) 

(.89. .95) 

I.88.J5) 

t .85. .03) 

1  too  hiah 

1  14(.o5i 

1 41.05) 

2  (..09 ) 

245(982) 

35(912) 

2o(.U'i 

i  cover 

!  2781.93) 

2531.84) 

2481.  S3) 

239(980) 

roo  low 

S(9»'l 

331.1 1) 

281.(10) 

l'(  .06) 

411.14) 

1  int.val 

j 

(.90. .96) 

(.S0..8S) 

(.'7, .86) 

(.'S..S') 

(."5. .84) 

too  hiah 

(Si  02' 

(S(.02) 

141.05) 

IO(  .06) 

10(.0' ' 

cover 

2S6<  .95) 

2Mi. 87) 

2'51.o2) 

I'-il.'M) 

2(0'. s') 

too  iow 

81.03  1 

331.11) 

1 1(904) 

'(9)2) 

'81  00) 

1  int.val 

,  (.93, .98) 

(•S3. .9 1 ) 

(.89.935) 

(.88, .95) 

(.84.91) 

TABLE  2 


LENGTH  OF  C.I  FOR  In  P(D>  t>  (EXPONENTIAL  MODEL) 
N  *  50,  0  =  0.5,  p/=  1.0,  p2=  10.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

SO  % 

0.2161 

(0.033) 

0.1902 

(0.032) 

0.0889 

(0.012) 

0.1546 

(0.080) 

0.1267 

(0.038) 

90  % 

0.2709 

(0.042) 

0.2447 

(0.041) 

0.1141 

(0.016) 

0.2048 

(0.106) 

0.1619 

(0.047) 

1.0 

SO  % 

0.3188 

(0.049) 

0.2882 

(0.042) 

0.1768 

(0.024) 

0.2022 

(0.066) 

0.1810 

(0.032) 

90  % 

0.4008 

(0.056) 

0.3717 

(0.055) 

0.2271 

(0.032) 

0.2680 

(0.088) 

0.2344 

(0.041) 

1.5 

SO  % 

0.4164 

(0.056) 

0.3798 

(0.053) 

0.2648 

(0.037) 

0.2822 

(0.088) 

0.2623 

(0.045) 

90  % 

0.5237 

(0.071) 

0.4919 

(0.079) 

0.3400 

(0.048) 

0.3740 

(0.117) 

0.3398 
(0.061)  • 

2.0 

80  0/o 

0.5238 

(0.082) 

0.4802 

(0.063) 

0.3528 

(0.049) 

0.3773 

(0.120) 

0.3522 

(0.065) 

90  % 

0.65S4 

(0.100) 

0.6253 

(0.076) 

0.4530 

(0.064) 

0.5001 

(0.159) 

0.4583 

(0.086) 

3.0 

SO  % 

0.7496 

(0.131) 

0.6925 

(0.125) 

0.5287 

(0.074) 

0.5846 

(0.188) 

0.5449 

(0.107) 

90  % 

0.9404 

(0.164) 

0.9173 

(0.178) 

0.6790 

(0.095) 

0.7749 

(0.249) 

0.7059 

(0.140) 

4.0 

80  % 

1.0498 

(0.260) 

0.9973 

(0.311) 

0.7047 

(0.099) 

0.7988 

(0.261) 

0.7380 

(0.148) 

90  % 

1.3136 

(0.323) 

1.4207 

(0.737) 

0.9049 

(0.127) 

1.0590 

(0.346) 

0.9589 

(0.194) 
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TABLE  3 

COVERAGE  RATIO  (EXP  MODEL) 
N  -  20,  0*0.5,  pj-  1.0,  p2-  10.0 


time  C.I  (%) 

coverage 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5  80  % 

(.7866) 

90  % 

too  high 
cover 
too  low 
int.val 

too  hish 
cover 
too  low 
int.val 

141.05' 

2721.91 

14(  .05 
(.S7..94) 

4(.01) 

282(.941 

14(.05 
(.9 1  ..97) 

401 .  3 
2281.76 
32.11 
(■7i;.81j 

401.13' 

2461.82) 

14(.05 

(-78..S6) 

241.08 
251  .84 

251.08' 
(.79, .88) 

14(.05' 

2761.92 

10(.03 

(.89,.95) 

106(.35) 

1741.58) 

201.07) 
(.52. .64) 

781.26) 

2111.70) 

,  11(041 
(.6a,.  /  6  > 

1401.47) 

160(.53) 

01.00) 

(.4S..59) 

991.33) 

20K.67) 

01.00) 

( .62.-72 ) 

1.0  SO  % 

(.6203) 

90  % 

too  hish 
cover 
too  low 
int.val 

too  high 
cover 
too  low 
int.val 

17(.06' 
250  .83* 

(.79..S8J 

2o(°92 

18(.06 

(.S8..95] 

401.13 
2271.76“ 
331.1  1 
(.71, .81) 

17(.06 
2651.88 
)  181.06 
(.85.. 921 

241 .08 
2501.83' 
261-09 
(.79, .88) 

141.05 
277(.92' 
91 .03  j 
(.89. .9a) 

721.24) 

2101.70) 

1S(.06> 

(.65,75) 

245(.S2) 

66(.  22  > 
2241.75) 
101.03) 

(.70, .80) 

331.1 1) 
259(.86) 
81.03) 

(.82. .90) 

1.5  80  % 

(.4891) 

90  % 

too  high 
cover 
too  low 
int.val 

too  hish 
cover 
too  low 
int.val 

151.05 

2651.88' 

20(.07 

(.S5..921 

15(.05 
278(.93 
71.02) 
(.88. .95) 

261.09 
226(.75 
481.16 
(.70. .80 

15(.05 
265(.8S 
201.07 
(.85;. 92' 

241.08)  481.16) 

250(.S3)  2471.82) 

26.09)  5(102) 

(.79,.88)  (.78..S7) 

14(.05)  331.11) 

2771.92)  2661.89) 

91.03)  11.00) 

(.S9..9a)  (.Sa..94) 

401 .13) 
2361.79)  ! 

24(  .08) 
(74..S3) 

191.06) 

2661.S9) 

151.05) 

(.85.-94)  | 

2.0  SO  % 

(.3857) 

90  % 

too  hish 
cover 
too  low 
int.val 

too  hish 
cover 
too  low 
int.val 

111.04 

w, 

(,85,.93j 

1 1  ( .  04 
2S01.93 

<« 

241.08 

2201.731 

56(.19 

(.68,.78j 

111.04 
2671.89 
22(.07 
(.85  .93 

24(  .08 
2501.83, 
26(.09 
(.79,.88j 

141.05' 

2771.92 

(.89«5;°93a 

M 

281.09) 

2711.90) 

(A®) 

271.09) 

2401.80) 

331.1 1) 
(.75,.85) 

101.03) 

2671.89) 

23  .08) 
(.S5..93) 

3.0  SO  % 

(.2399) 

90  % 

too  hish 
cover 
too  low 
int.val 

too  hish 
cover 
too  low 
int.val 

231.08)  231.08 

271(.90)  2491.83 

61.02)  28(.09 

(.S7..94)  (.79,.87] 

81.03)  8(.03) 

286(.9a)  2641.88 

61.02)  28  .09 

(.93. .98)  1.84. .921 

231 .08 
25  l(.S4 
26.09, 
(.79,. 88) 

14(.05' 

2771.92 

91.03J 

(.S9..9a 

401 .13) 
252  .84) 
81.03) 
(.80,  .88) 

22(.07) 

2761.92) 

l.S^^/a) 

191.06) 

238  79) 

43  .14) 
(75,-84) 

M<\ 

1.83..91) 

4.0  80  % 

(.1492) 

90  % 

1 _ 

too  hish 
cover 
too  low 
int.val 

too  hish 
cover 
too  low 
int.val 

191.06, 

2711.90 

101.03 
(.87, .94) 

31.01) 

2871.96 

W.D2) 

(.93..9S) 

191.06 

2241.75) 

57(.19 
(.70, .80) 

31.01) 

2401.80 

281.09 

(.75..S5) 

23<.OS 
251(.84' 
26.09 
(.79,  .88) 

14(.05 
27  (.92 
9)  .03) 

(.89,-95. 

253^S4) 
111.04) 
(.8 1..89) 

211.07) 

2  1.  >2) 
:i.on 
(.S9..95) 

151.05) 

231  .77 
541,18) 

(72. .82) 

71.03) 

36  f(  .S*  i 
321.11) 

(.83. .91)  | 

TABLE  4 

LENGTH  OF  C.I  FOR  In  P(D<t}  (EXPONENTIAL  MODEL) 
N  -  20,  6 -  0.5,  Pj  *  1 .0,  p, *  10.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.3696 

(0.092) 

0.3026 

(0.092) 

0.1446 

(0.033) 

0.2083 

(0.079) 

0.1854 

(0.059) 

90  % 

0.4598 

(0.115) 

0.3897 

(0.122) 

0.1858 

(0.042) 

0.2713 

(0.103) 

0.2449 

(0.081) 

|  1.0 

80  % 

0.5451 

(0.131) 

0.4703 

(0.125) 

0.2877 

(0.066) 

0.3072 

(0.108) 

0.2899 

(0.081) 

90  % 

0.6797 

(0.164) 

0.6132 

(0.168) 

0.3695 

(0.085) 

0.4001 

(0.140) 

0.3816 

(0.110) 

1.5 

80  % 

0.7156 

(0.167) 

0.6320 

(0.159) 

0.4308 

(0.099) 

0.4668 
(0.1  SO) 

0.4361 

(0.130) 

I 

i 

90  % 

0.8928 

(0.208) 

0.8351 

(0.227) 

0.5533 

(0.127) 

0.6079 

(0.235) 

0.5790 

(0.181) 

I  2.0 

SO  % 

0.9150 

(0.273) 

0.8328 

(0.306) 

0.5739 

(0.132) 

0.6452 

(0.25S) 

0.5971 

(0.186) 

90  % 

1.1412 

(.0.340) 

1.1847 

(0.837) 

0.7370 

(0.169) 

0.8403 

(0.336) 

0.7927 

(0.259) 

3.0 

SO  % 

1.3762 

(0.567) 

1.3689 

(0.878) 

0.8601 

(0.198) 

1.0177 

(0.419) 

0.9324 

(0.307) 

90  % 

1.7105 

(0.663) 

2.2324 

(1.773) 

1.1045 

(0.255) 

1.3253 

(0.546) 

1.2427 

(0.429) 

4.0 

80  % 

2.1097 

(1.197) 

2.4656 

(2.060) 

1.1463 

(0.264) 

1.3966 

(0.584) 

1.2722 

(0.426) 

1 

90  % 

2.5832 

(1.308) 

3.9914 

(2.720) 

1.4720 

(0.339) 

1.8188 

(0.761) 

1.6956 

(0.596) 
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TABLE  6 

LENGTH  OF  C.I  FOR  In  P(D>  t}  i EXPONENTIAL  MODEL) 
N  -  50,  6-2/3,  pj -1.0,  p2-  10.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.2536 

(0.040) 

0.2261 

(0.038) 

0.1196 

(0.017) 

0.2232 

(0.092) 

0.1914 

(0.048) 

90  % 

0.3184 

(0.050) 

0.2913 

(0.049) 

0.1535 

(0.022) 

0.2958 

(0.122) 

0.2506 

(0.060) 

1.0 

80  % 

0.3848 

(0.056) 

0.3502 

(0.052) 

0.2383 

(0.034) 

0.2618 

(0.082) 

0.2510 

(0.048) 

90  % 

0.4840 

(0.070) 

0.4530 

(0.069) 

0.3060 

(0.044) 

0.3470 
(0. 109) 

0.3243 

(0.064) 

1.5 

80  % 

0.5260 

(0.085) 

0.4823 

(0.080) 

0.3570 

(0.051) 

0.3841 

(0.134) 

0.3635 

(0.074) 

9<  % 

0.6612 

(0.107) 

0.6282 

(0.107) 

0.4584 

(0.066) 

0.5091 

(0.178) 

0.4783 

(0.102) 

2.0 

80  % 

0.6815 

(0.126) 

0.6284 

(0.120) 

0.4757 

(0.068) 

0.5387 

(0.196) 

0.4977 

(0.106) 

90  % 

0.8556 

(0.157) 

0.8281 

(0.173) 

0.6108 

(0.088) 

0.7140 

(0.260) 

0.6584 

(0.144) 

3.0 

80  % 

1.0966 

(0.332) 

1.0598 

(0.468) 

0.7131 

(0.102) 

0.8739 

(0.332) 

0.7S67 

(0.175) 

1 

1 

90  % 

1.3718 

(0.414) 

1.5369 

(0.930) 

0.9157 

(0.132) 

1.1583 

(0.440) 

1.0428 

(0.237) 

4.0 

80  % 

1.7374 

(0.773) 

1.9119 

(1.277) 

0.9505 

(0.137) 

1.2188 

(0.476) 

1.0833 

(0.248) 

i 

90  % 

2.1575 

(0.901) 

3.1281 

(1.968) 

1.2205 

(0.176) 

1.6155 

(0.631) 

1.4384 

(0.336) 

i 


i 


time  C.I  (%)  coverage 


0.5  80  %  too  high 

cover 
too  low 
int.val 


NOR 


MLE 


i("231) 


90  % 


too  high 
cover 
too  low 
int.val 


(.5241) 


)  31l 

[.io: 

1  33 

i  23  li 

.77] 

)  233 

i  3  Si 

1.131 

l  34 

(.70. 

,82) 

(-73 

l  32l 

r.  ii' 

)  18 

►  253i 

.841 

i  268 

i  161.05 

l  14 

(.SO, 

..88) 

(.86 

(.3799) 


too  man 
cover 
too  low 
int.val 

too  high 
cover 
too  low 
int.val 


32f.ll) 
232C.77 ) 
36(.  12) 
(.73,-82) 

I4f.05) 
273  .91) 
13  .04) 
(.S8..941 


21(.07) 
251  ,S4 
28(.09) 
(.79,88) 

7C.02) 

2871.96) 


32f.ll) 


141.05) 
250(.S3) 
36( .  1 2) 
(.79..S8) 


211.07) 
2211.74) 
58(.  19) 
(.69.79) 


2651. 88) 
28(.09) 
(.85.92) 


j  (.  1446 ) 


4.0  8 


(.0760) 


321.11) 
2341. -S) 
34(.  1 1) 
(.73,83) 

1 8(  .06) 

2681.89) 
141.05) 
(.86,93) 


321.1 1 ) 
235  .78) 
33(.ll) 
(.74,. 83) 

18(06) 

2691.90) 

13(.04) 

(.86.93) 


321.11) 
2351.78) 
33(.l  1) 
(.74,. 83) 

18(.06) 

2701.90) 
12  .04) 
(.S7..93) 


101(.34) 

1931.64) 

61.02) 

(.59..70) 

691.23) 

2301.77) 

11.00) 

(.72,81) 


811.2^) 

2111.70) 

81 .03 ) 
(.65, .76) 

52(.17) 

2451.82) 

31.01) 
(.77,. 86) 


64(.21) 
227  .76) 
91.03) 
(.71, .81) 

33(.l  1) 
2621.87) 
51.02) 
(.83, .91) 


501.17) 
239(.80) 
11  .04) 
(.75,-84) 

28(.09) 
2671.89) 
5(.02) 
(.So, .93) 


B.T 


222(.74 

77(.26 

11.00) 
(.21. .31 

1871.62 
1131.38 
01.00) 
(.32, .43' 


1121.37) 

176(.59) 

12(.04) 
(.53, .64) 

761.25) 

2171.72) 

71.02) 
(.67. .77) 


671.22) 


261.09) 

64,74) 

361.12) 

241(.80) 

23(.08) 
(.76..  35) 


38(.13) 
227  .76) 


227  .76) 
35(.12) 
(.71,81) 

15(.05) 

2631.SS) 

221.07) 

(.84,92) 


161.05 
237(.79‘ 
47  .16' 
(.74,84; 

281.09' 

(.S4.‘.93l 


15(.05 
2311.77 
54(.  IS 
(.71,82 

61.02) 
2611.87 
33(.  1 1 
(.82,92 


TABLE  8 


LENGTH  OF  C.I  FOR  In  P{D>t}  (EXPONENTIAL  MODEL) 
N  =  20,  0»2  3,  px  —  1.0,  p2-  10.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

80  % 

0.4304 

(0.109) 

0.3617 

(0.104) 

0.1981 

(0.053) 

0.3451 

(0.172) 

0.3168 

(0.128) 

90  % 

0.5359 

(0.173) 

0.4680 

(0.139) 

0.2544 

(0.068) 

0.4494 

(0.225) 

0.4302 

(0.183) 

1.0 

80  % 

0.6594 

(0.170) 

0.5788 

(0.161) 

0.3952 

(0.106) 

0.4630 

(0.207) 

0.4394 

(0.164) 

90  % 

0.8227 

(0.212) 

0.7618 

(0.225) 

0.5075 

(0.136) 

0.6030 

(0.298) 

0.5937 

(0.242) 

1.5 

80  % 

0.9426 

(0.323) 

0.8710 

(0.469) 

0.5931 

(0.161) 

0.7228 

(0.399) 

0.6560 

(0.256) 

90  % 

1.1756 

(0.402) 

1.2131 

(0.773) 

0.7616 

(0.207) 

0.9413 

(0.519.) 

0.8865 

(0.379) 

2.0 

so  % 

1.2684 

(0.587) 

1.2540 

(0.S99) 

0.8300 

(0.222) 

1.0229 

(0.603) 

0.9008 

(0.359) 

90  % 

1.5765 

(0.68S) 

2.0258 

(1.817) 

0.9311 

(0.2S6) 

1.3322 

(0.785) 

1.2240 

(0.517) 

3.0 

| 

SO  % 

2.2188 

(1.311) 

2.5096 

(2.079) 

1.1951 

(0.400) 

1.6583 

(1.021) 

1.4301 

(0.574) 

1 

90  % 

2.6995 

(1.400) 

4.1815 

(2.807) 

1.5347 

(0.914) 

2.1597 

(1.330) 

1.9457 

(0.830) 

80  % 

3.5128 

(1.957) 

3.1967 

(2.724) 

1.6119 

(0.778) 

2.3080 

(1.443) 

1.9704 

(0.796) 

1 

90  % 

4.0702 

(1.893) 

4.8194 

(3.183) 

2.0699 

(0.993) 

3.0057 

(1.SS0) 

2.6794 

(1.140) 

C.  ROBUSTNESS 

The  robustness  of  the  estimates  in  Chapter  II  was  studied  with  respect  to  an 

incorrect  model  assumption  about  the  distribution  of  the  sojourn  time  in  state  1.  In 

the  previous  simulations,  the  maximum  likelihood  estimator  used  the  known  correct 
model. 

The  data  for  the  simulation  experiment  in  this  section  is  generated  from  the 
following  three  state  semi-Markov  process:  Individuals  start  in  state  1  at  t  =  0.  The 
probability  of  a  transition  to  state  0  is  0,  and  to  state  2  is  1-0.  From  state  2,  the 

probability  of  a  transition  to  state  l  is  1.  State  0  is  an  absorbing  state.  The  sojourn 


time  in  state  2  is  exponential  with  mean  l/p2.  The  sojourn  time  in  state  1  is  the  sum  of 
two  independent  exponentials  with  means  l/p1  and  l/p3;  that  is,  the  sojourn  time  in 
state  1  has  a  hvpoexponential  distribution.  The  same  Fortran  program  for  the 
simulation  is  used  but  slightly  modified  for  the  above  change.  The  data  generated  are 
analyzed  by  the  same  Fortran  programs  for  each  estimator  as  in  the  first  section.  In 
particular,  the  maximum  likelihood  estimator  assumes  the  sojourn  time  in  state  1  has 
an  exponential  distribution  rather  than  the  true  hypoexponentiai  distribution. 

For  the  first  simulation  results  reported  in  Tables  9  to  12.  parameter  values  of  pt 
=  2.0,  p,  =  10.0,  p3  =  2.0,  and  0  =  0.5  are  used.  Again  two  different  numbers  of 
individuals  are  used;  they  are  20  and  50.  The  simulation  is  replicated  300  times,  the 
coverage  numbers  and  the  average  lengths  of  the  confidence  intervals  are  computed. 
The  actual  value  of  the  survival  function  is  computed  by  inverting  the  Laplace 
transform  of  the  passage  time  to  state  0  for  the  semi-Vlarkov  process.  Some  computed 
values  of  the  survivor  function  can  be  found  in  Appendix  B. 

Table  9  reports  coverage  results  for  the  case  N  =  50  individuals.  The  use  of  the 
incorrect  model  for  the  maximum  likelihood  estimator  leads  to  the  majority  of  the 
intervals  being  too  low  for  t  =  0.5  and  t*  1.0.  The  binomial  confidence  interval  has  a 
slight  tendency  to  overcover.  The  confidence  interval  based  on  the  asymptotic  renewal 
estimator  undercover  for  t  —  0.5  and  t*  1.0.  For  larger  values  of  t  they  have  the 
correct  coverage.  The  results  for  the  jackknife  and  bootstrap  individuals  are  similar  to 
those  in  Table  1.  Table  10  reports  the  average  lengths  of  the  confidence  intervals.  The 
average  lengths  for  the  maximum  likelihood  procedure  are  similar  to  those  of  Table  2. 
For  the  other  procedures,  the  average  lengths  are  smaller. 

Results  of  the  simulation  expenment  for  N  =  20  individuals  are  given  in  Table 
11  and  12.  Table  11  shows  the  coverage  results.  The  maximum  likelihood  procedure 
undercovers  for  t  =  0.5  and  t=  1.0  and  overcovers  for  t=  1.5  and  t  =  2.0.  The  binomial 
procedure  once  again  tends  to  overcover.  The  procedures  based  on  the  asymptotic 
renewal  estimator  have  the  correct  coverage  for  all  t  except  t  =  0.5.  The  average  lengths 
of  the  intervals  are  given  in  Table  12. 

Tables  13  to  16  report  results  of  simulation  of  a  semi-Markov  model  with  0=  2:3 
but  with  the  other  parameters  the  same;  p;  =  2.0,  p,=  10.0,  p,=  2.0  . 

Results  for  N  =  20  individuals  appear  in  Tables  15  and  16.  The  confidence 
intervals  based  on  the  asymptotic  renewal  estimator  have  the  correct  coverage  for  all  t 
except  0.5  and  1.0.  The  resuits  are  similar  to  those  of  Tables  11  and  12. 


TABLE  9 


CO VERA 
N  -  50, 


RAGE  RATIO  (HYPOEXP  MODEL) 
50,  0-0.5,  pj  -  2.0,  p2=  10.0,  p3  =  2.0 


time  C.I  (%) 


0.5  SO  % 


(.8652) 


length 

BIN 

NOR 

MLE 

J.K 

B.T 

too  high 
cover 
too  low 
int.val 

1S6.06) 

261(.S7) 

21  (.07) 
(.S3..91) 

244|:sf) 

21  (.07) 
(.77..S6) 

06.00) 

36.01) 

297(.99) 

(.00..01) 

130(.43) 

1656.55) 

56.02  V 
(.49..01) 

1496.50) 

1506.50) 
l(.OO) 

!  .44. .56) 

90  % 


too  high 
cover 
too  low 
int.val 


106.03) 

283(.94) 


1 8(  .06) 
275(.92) 


0(.00) 

6(,02) 

294(.9S) 
( .00. .04) 


76.02J  76.02) 

.92. .97)  (.S9..90) 


99(.33) 

2006.67) 


106^  .35 ) 
1946.65) 
O(.OO) 
(.59. .70) 


1.0  SO  % 

(.6743) 

too  high 
cover 
too  low 
int.val 

226.07) 

2606.87) 

186.06) 

(.S3..91) 

36.01) 
2296.76) 
336.1  1) 
(-72..S1) 

26.01) 

145(.48) 

153(.5D 
(.43, .54) 

426.14 
242(.8T 
166.05 
(.76, .85' 

396.13) 

2406.  SO) 

216.07) 

(.75,  .85) 

90  % 

too  high 
cover 

116.04) 

280(.93) 

226.07) 

2696.90) 

06.00) 

:&:31 

6.662.76) 

286.09 

2636.S8 

216.07) 

2656.88) 

too  low 
int.val 

96.03) 

6. 91. .96) 

.96.03) 
(.86. .93) 

9(  .03) 

6 . 84  ..91) 

146.05)  -i 

6  .S5..92) 

(.5157) 


too  high 
cover 
too  Low 
int.val 


276.09 

254(.S5 

196.06 


)  276.09) 

66.02) 

316.10) 

256.08) 

2356.78) 

2426.81) 

247(.S2) 

2436.81) 

)  3S(.  13) 

526.17) 

226.07) 

32(.ll) 

)  (.74.83) 

(.76..851 

(.7S..S7) 

(.77..S5) 

too  high 
cover 
too  low 
int.val 


1 1(.04) 
2766.92) 


2766.92) 
13  (.04) 
(.89, .95) 


21(.07) 
261  (.87) 
18  .06) 
6.S3..91) 


4(.01) 

276(.93 

201.07 


201 .07 
(.S9..95 


216.07) 
271  (.90) 
S6.02) 
(.87. .94) 


19(.06) 
62  .87) 


262  .87) 
19  .06) 
(.84..91) 


!  2.0  8 


too  high 
cover 
too  low 


(.3930) 


256.08 
536.84 
22  .07 
80.. 88 


25(.08 

240(.S0 

35(.12 

(.75.85' 


too  high 
cover 
too  low 
int.val 


166.05) 
2746.91) 
106.03) 
(.88..95) 


166.05) 
2626.87) 
22(.07) 
(.84..91 ) 


6.02) 

282(.94) 

116.04) 
( .9 1  ..97) 


20(.O7) 
2696.90) 
1 1  ( .04) 
6. 86. .93) 


126.041 
2646.8  S) 
246.081 
(.84. .92) 


(.2278) 


too  high 
cover 
too  low 
int.val 


25(.08) 

2476.82) 


(.78..S7 


256.08) 

220(,73) 

55(.1S) 
(.68, .78) 


346.1  1) 
2536.84) 
13(.04) 
(.80, .88) 


316.10) 

2486.83) 

216.07) 

(.7S..S7) 


216.07) 
2416.80) 
386.13) 
(."6. .85 1 


90  %  |  too  high 

176.06) 

17(.06) 

20(.07) 

1S(.06) 

lt( 

1  cover 

273  .91) 

2556.85) 

2756.92) 

27 1  ( .90) 

260 

!  too  low 

10  .03) 

/  t  >  n  i  •»  \ 

2S(.09) 

56.02) 

i  Of)  .  1  ?  \ 

11  .04) 

/  at'  n  ii 

236 

(.89. .95)  ( .87. .94) 


too  high 
cover 
too  low 
int.val 

116.04) 

2666.89) 

236.08) 

(.85,.92) 

11  (.04) 
235  .78 
546.18) 
(.74, .83) 

47(.  1 6 1 
246  .82) 
"6.02) 
(.78, .86) 

3K.10) 
2446.81) 
256.08) 
(.77, .86)) 

206.07) 
23S6."9) 
4  26.14) 
( .7  >..84) 

too  high 
cover 
too  low 
int.vai 

66.02) 
2856.95) 
96.03) 
(-93,  .97) 

66.02) 

27(6.90) 

236.08) 

(.87„94) 

296.10) 

2676.89) 
41.01  i 
(.85, .93) 

196.06) 
2~o(  oo) 

1 16.04) 
(.87, .93) 

106.03) 
2  ('46. 88) 
2(m.oo) 
(.84, .92) 

ww 


TABLE  10 

LENGTH  OF  C.I  FOR  In  PJD  >  t}(HYPOEXP  MODEL! 
N  -  50.6  *  0.5.p  j  -  2.0, p2  =  10.0, p3  =  2.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

b:t 

0.5 

80  % 

0.1688 

(0.029) 

0.1442 

(0.029) 

0.0886 

(0.011) 

0.0913 

(0.042) 

0.0806 

(0.022) 

90  % 

0.2110 

(0.037) 

0.1852 

(0.037) 

0.1138 

(0.014) 

0.1210 

(0.056) 

0.1043 

(0.028) 

1.0 

80  % 

0.2843 

(0.043) 

0.2553 

(0.041) 

0.1764 

(0.022) 

0.1666 

(0.053) 

0.1534 

(0.028) 

90  % 

0.3572 

(0.054) 

0.3291 

(0.053) 

0.2265 

(0.028) 

0.2208 

(0.071) 

0.1988 

(0.037) 

1.5 

80  % 

0.3958 

(0.058) 

0.3606 

(0.054) 

0.2641 

(0.033) 

0.2692 

(0.083) 

0.2496 

(0.047) 

90  % 

0.4978 

(0.0731 

0.4666 

(0.072) 

0.3392 

(0.042) 

0.3568 
(0. 110) 

0.3240 

(0.064) 

2.0 

so  % 

0.5139 

(0.080) 

0.4710 

(0.075) 

0.3519 

(0.044) 

0.3777 

(0.117) 

0.3506 

(0.067) 

90  % 

0.6461 

(0.100) 

0.6130 

(0.100) 

0.4519 

(0.056) 

0.5006 

(0.155) 

0.4561 

(0.093) 

3.0 

80  % 

0.7840 

(0.160) 

0.7261 

(0.157) 

0.5274 

(0.066) 

0.5998 

(0.188) 

0.5576 

(0.111) 

90  % 

0.9832 

(0.199) 

0.9684 

(0.247) 

0.6773 

(0.085) 

0.7950 

(0.249) 

0.7261 

(0.153) 

4.0 

80  % 

1.1592 
(0  334) 

1.1276 

(0.480) 

0.7030 

(1.088) 

0.8242 

(0.260) 

0.7650 

(0.154) 

90  % 

1.4496 

(0.415) 

1.6667 

(1.017) 

0.9027 

(1.113) 

1.0924 

(0.345) 

0.9980 

(0.214) 

TABLE  11 


COVERAGE  RATIO  (HYPOEXP  MODEL) 
N  -  20.  6 -  0.5,  Pj  -  2.0,  p,  -  10.0,  p^  =  2.0 


P21 


Py 


time  C.I  (%) 


0.5  80  % 


(.$652) 


1.0  80  % 


(.6743) 


90  % 


1.5  80  0/o 


(.5157) 


90  % 


2.0  80  % 


(.3930) 


90  % 


3.0  SO  % 


(.2278) 


90  % 


4.0  80  % 


(.1320) 


length 


too  high 
cover 
too  low 
int.val 

too  high 
cover 
too  low 
int.val 


too  high 
cover 
too  low 
int.val 

too  high 
cover 
.too  low 
int.val 


too  high 
cover 
too  low 
int.val 

too  high 
cover 
too  low 
int.val 


too  high 
cover 
too  low 
int.val 

too  high 
cover 
too  low 
int.val 


too  high 
cover 
too  low 
int.val 

too  high 
cover 
too  low 
int.val 


90  % 


too  high 
cover 
too  low 
int.val 

too  high 


BIN 

NOR 

MLE 

J.K 

B.T 

12(.04) 
2781.9  3) 

10(.03) 
(.90, .96) 

591.20' 
231  (.IT 
10.03' 
(.72,. 821 

1  01.00) 

>  531.18) 

)  247(.S2) 

1  (.13, .22) 

1051.35) 

1831.61 

12(.04 
(.55, .67) 

i— 

iL. 

O _ 

01.00) 
290(,97) 
10(.03) 
(.95  .99) 

:Si- 

(.83;092) 

)  01.00) 

)  109(.3o) 

1 9 1  ( .  64) 

'  (.31.42) 

>  (.6ll;°3 % 

161.05) 

254(.S5) 

30(.10) 

(.80,.88) 

51.02) 

2S5(.95) 

10(.03) 

(.93.97) 


44(.15) 
2261.75) 
30  .10) 
(.70,  .80) 

16(.05) 
274  .91 
101.03) 
(. 881.95) 


11.00) 

2211.74) 

78(.26) 

(.69,.79) 

01.00) 
266(.S9) 
34  .11) 
(.85.. 92) 


45(.  15) 
2401.80 
15(.05) 
(.7  5..S5) 

291.10) 

263(.88) 


321.11) 
241  .80) 
27  .09) 
(.76, .85) 

1 1  (.04) 
271  .90) 
181.06) 
(.S7..94) 


211.07) 

260(.S7) 

19(.06) 

(-S3..91) 

51.02) 

2761.92) 

191.06) 

(.S9..95) 


211.07) 
2421.S1) 
37  .12) 
(.76..S5) 

21(.07) 

2601.87) 

191.06) 

(.83;.9t) 


81.03) 
257(.86) 
35(.  12) 
(.82..90) 

21.01) 

2831.94) 

151.05) 

(.92.97) 


401.13) 
2441.81) 
16(.05l 
(.77,  .86) 

24(.08) 
273  .91 
31.01) 
(.S8..94) 


231 ,0S) 
2321.77) 
45(.15) 
(.3, .82) 

7(.02) 

264(.S8) 

291.10) 

l.84;.92) 


1 1(.04) 
2701.90) 
19(.06) 
(.87,.93/ 

1 11 .04) 
2841.95) 
o(.02) 
(.91.97) 


28(.09) 
235  .78) 
37  .12) 
(.74,.83) 

1 1(04) 
2701.90) 
19  .06) 
(.87.93) 


171.06) 

2601.87) 

23(.OS) 

(.83,-91) 

71.02) 

284(.95) 

91.03) 

(.92..9/) 


391.13) 

2461.82) 

15(.05) 

(.78,-86) 

191.06) 
277  .92) 
41.0  lj 
(.89,  .95) 


18(06) 
232  .77) 
50(.17) 
(.73,.82) 

263i°8'8) 
30  .10) 
(.84, .91) 


131.04) 

2731.91) 

14(.05) 
(.88, .94) 


51.02) 

281(.94 

141.05 
1.91. .96 


13(.04) 
2491.83) 
3S  .13) 
(.79,.87) 

51.02) 
257(,86) 
38(.  13) 
(.82.. 90) 


311.10) 
2541.85) 
15  .05) 
(.81,-89) 

17(.06) 
279  .93) 
4(.0h 
(.90, .96) 


37(.  12) 
251  .S4) 
12  .04) 
(.79.88) 


171.06) 
2301.77 
53(.18 
(• 72..81 

61.02) 
26 1( .  S7) 
33  .11) 
1.83. .91) 


41.01) 
279(.93) 
171.06) 
(.87, .94) 

41.01) 


22^. ll) 

701.23) 
(.70,.  80) 

21.01 ) 


44(.15) 
245  ,S2) 
II  (.04) 
(.79, .88) 

231.08) 


3S(.  13) 
2511.S4) 
111.04) 
(-81,-89) 

191.06) 


1 7(  .06 ) 
229  .76) 
541.18) 
1-72,  .82) 


TABLE  12 

LENGTH  OF  C.I  FOR  In  P(D>  t!(HYPOEXP  MODEL) 
N  -  20,6  -  0.5, p  j  -  10.p2  -  10.0, p3  -  2.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  0/o 

0.2891 

(0.077) 

0.2217 

(0.083) 

0.1438 

(0.030) 

0.1405 

(0.056) 

0.1296 

(0.046) 

90  0/o 

0.3591 

(0.095) 

0.2810 

(0.111) 

0.1846 

(0.03S) 

0.1829 

(0.075) 

0.1750 

(0.066) 

1.0 

80  % 

0.4845 

(0.113) 

0.4131 

(0.107) 

0.2861 

(0.059) 

0.2694 

(0.105) 

0  *“07 
(0.073) 

90  0/o 

0.(5038 

(0.142) 

0.5365 

(0.143) 

0. 3674 
(0.076) 

0.3509 

(0.136) 

0.3314 

(0.102) 

1.5 

80  % 

0.6810 

(0.175) 

0.5995 

(0.168) 

0.4284 

(0.088) 

0.4512 

(0.180) 

0.4190 

(0.12S) 

90% 

0.8496 

(0.218) 

0.7911 

(0.241) 

0.5501 
(0;  113) 

0.5876 

(0.234) 

0.5517 

(0.174) 

o 

r  i 

80  % 

0.8843 

(0.230) 

0.7974 

(0.240) 

0.5707 

(0.118) 

0.6432 

(0.259) 

0.5944 

(0.290) 

90  % 

1.1030 

(0.286) 

1.0914 

(0.513) 

0.7329 

(0.151) 

0.8376 

(0.337) 

0.7864 

(0.256) 

3.0 

80  % 

1.4552 

(0.539) 

1.5826 

(1.296) 

0.S553 

(0.176) 

1.0345 

(0.419) 

0.9528 

(0.313) 

90  % 

1,8145 

(0.737) 

2.48S1 
( I.9S4) 

1.0983 

(0.227) 

1.3472 

(0.546) 

1.2605 

(0.427) 

4.0 

so  % 

2.3431 

(1.375) 

2.6289 

(2.141) 

1.1399 

(0.235) 

1.4289 

(0.581) 

1.3138 

(0.437) 

90  % 

2.8401 
i  1.447) 

4.3272 

(2.778) 

1.4638 

(0.302) 

1.8610 

(0.757) 

1.7390 
( 0.60 1 ) 

TABLE  13 


COVERAGE  RATIO  (HYPOEXP  MODEL) 
N  =  50,0-2:3,  pj  =2.0,  p2-  10.0, p,  -  2.0 


time  C.I  (%) 

length 

BIN 

NOR 

MLE 

J.K 

B.T 

too  hieh 

25(.08) 

251. OS) 

01.00) 

1S61.62) 

221.1.73) 

cover 

2491.83) 

2491.83) 

21.01) 

1091.36) 

801.27) 

too  low 

26(.09) 

261.09) 

2981.99) 

51.02) 

O1.0O) 

int.val 

(.78. .87) 

(.78. .87) 

(.00. .02) 

(.3 1..42) 

( .22. .32) 

(.8214) 

too  hieh 

81.03) 

251 .08) 

01.00) 

41.01) 

2961.99) 

1571.52) 

1941.65) 

cover 

284(.95) 

2671.89) 

1411.47) 

1061.35) 

too  low 

81.03) 

SI. 03) 

21.01) 

01.0(0 

int.val 

(.92. .97) 

(.85. .93) 

(.00. .03) 

(.41. .53) 

(.30. .4 1 ) 

1.0  so  % 

too  hieh 

301.10) 

301.10) 

01.00) 

651.22) 

591.20) 

cover 

2541.85 

2441.81) 

1791 .60) 

219(.73) 

225(.76) 

too  low 

16(.05) 

261.09) 

1211.40) 

11  (.04) 

lo<  .05) 

int.val 

(.SI. .89) 

(.77..S6) 

1.54,.o5) 

(.6S..7S) 

(.71. .81) 

(.5788) 

90  % 

too  hieh 

1 0(  .03) 

151.05) 

01.00) 

391.13) 

321.11) 

cover 

2801.93) 

2691.90) 

237(.79) 

2521.84) 

2611.87) 

too  low 

101.03) 

161.05) 

631.21) 

.91.03) 

’(.02) 

int.vai 

( .9 1  ..97) 

(.S6..93) 

1. 74.284) 

(.SO. .88) 

1.S3. .91) 

(.3936) 


i  9Q  o', 


2.0  80  % 


!  (.2652) 


90  % 


too  hieh 
cover 
too  low 
int.val 

191.06) 

2531.84) 

281.09) 

(.802.88) 

191.06) 

2421.81) 

391.13) 

(.761.85) 

81.03) 

2641.88) 

281.09) 

(.842.92) 

421.14) 

2411.80) 

171.06) 

(.762.85) 

291.10) 
2401.80) 

311.10) 
(.752.85) 

too  hieh 
cover 
too  low 
int.val  ! 

111.04) 
2741.9 1 ) 
151.05) 
1.S8..95) 

111.04) 
2611.87) 
281.09) 
(.83. .91) 

41.01) 

2841.95) 

121.04) 

(.922.97) 

171.06) 
2751.92) 
81.03) 
1.89. .95) 

151.05) 
265i.S8) 
20(.07) 
(.85..  92) 

too  hieh 
cover 
too  low 
int.val 


too  hieh 
cover 
too  low 
int.val 


261.09) 

2491.83) 

25(.U8) 

(.79,.87i 


51.02) 

283(.94) 

121.04) 

(.92..97) 


141.05) 
261  (.87) 
25(.08) 
(.83, .9 1 ) 


11  (.04) 
285(.95) 


.95) 
4(.on 
(.93.  .9') 


18(.06) 
2T4(  .9 1 ) 
S(.03) 
(.88, .95* 


I2(.04) 
262(.87) 
26(.<)9) 
(.84, .91 ) 


>1 .09) 

271.09) 

351.12) 

241.08) 

(.78) 

26  If.  87) 

249(.83) 

2351. 78) 

(.14) 

12(.04) 

16  .07) 

4 1  ( . 1 4) 

'.82) 

(.83,-91) 

(.79.87) 

(.74.. S3) 

i  3.0  SO  'Vo 

(.1197) 

90  % 

too  hieh 
cover 
too  low 
int.val 

j  too  hieh 
cover 
too  low 
int.val 

141.05) 

271(.90) 

15(.U5) 

(.872.94) 

51.02) 
2921.97) 
31.01) 
(.96. .99) 

141.05) 
2431.81) 
43(.  14) 
(.77,.85) 

51.02) 

1.802.88) 

691.23) 

2271.76) 

41 20 1 ) 
(.71, .81) 

37(.12) 
2611.87) 
21.01) 
(.81.9  0 

351.12) 
2491.83) 
16(  .00) 
(.79, .87) 

1 7(.06) 
2751.92) 
8(203) 
1.S6..95) 

221.07) 
2291.76) 
491.16)  j 

(.72. .81)  j 

61.02J  j 

2611.8  ) 
331.11) 

(.832.90 

4.0  80  % 

too  hieh 

131.04) 

131.04) 

981.33) 

341.11) 

191.06)  i 

cover 

275(.92) 

2211.74) 

1991.66) 

24S» 

r.83) 

2291.76) 

too  low 

121.04) 

661.22) 

31.01 ) 

ISl 

i  .06) 

521.1") 

(.0540) 

int.val 

1  (.89. .95) 

(.69..  79) 

(.61. .72) 

(.78, 

..87) 

I.-2..SD  i 

i 

9()  % 

too  hieh 

!  41.01) 

41.01) 

631.21 ) 

151.05) 

61.02)  1 

! 

cover 

2961.98) 

2301.-7) 

2371.-9) 

2— i 

’.92) 

Z'y  .s  " ) 

too  low 

0(  .1)0) 

661.22) 

D(  OO) 

si. 

O’) 

3-il.  10 

int.val 

(.97, .99) 

(.72. .81) 

(.74. .84) 

(.89, 

,.95) 

1.S3..91)  j 

TABLE  15 


COVERAGE  RATIO  (HYPOEXP  MODEL) 
N  =  20,6  =  2/3,  pj  - 2.0,  p2=  10.0, p3  =  2.0 


time  C.I  (%) 


0.5 


(.8214) 

'  A 


70 


1.0  80  % 


(.5788) 


90  % 


1.5  80  % 


j  (.3936) 


90  % 


3.0  80  % 


(.1197) 


90  % 


4.0  80  % 


(.0540) 


90  % 


length 


too  high 
cover 
too  low 
int.val 

too  high 
cover 
too  low 
int.val 


too  high 
cover 
too  low 
int.val 

too  high 
cover 
too  low 
int.val 


2.0  80  % 

too  high 
cover 
too  low 
int.val 

(.2652) 

too  high 

90  % 

cover 

too  low 
int.val 

BIN 

NOR 

MLE 

J.K 

B.T 

31.01) 

277(.92) 

20(.07) 
(.89, .95) 

261.09) 
254  .S5) 
20  .07) 
(.81, .88) 

01.00) 

541.18] 

246(.82 

(.14,.22) 

1241.41) 

1  172(.57) 

)  41.01) 

i  (.52, .63) 

1981.66) 

101(-34) 

11.00) 

(.2S..39) 

28&?96) 

91.03) 

(-94..9S) 

26(.09) 

265(.8S) 

(.si:°93i) 

01.00) 
1  lo(.381 
1 85<.62 
(.33.441 

881.29) 

1  2111.70) 

!  (.6x7^) 

1441.48) 

156(.52) 

01. 00) 
(-28,-39) 

25(.08) 
242  .81 
33(.l  1) 
(.76, .85) 

131.04) 

267(.S9) 

201.07) 

(.85..931 


251.08) 
242  .81) 
33(.l  1) 
(.76, .85) 

221.07) 
2601. 87) 
18  .06) 
(.83  .91) 


01.00) 
238(.79 
62(.21 
(.75, .84 


0(.00) 
2741.91) 
26  .09) 
(.88,.95) 


371.12) 

2531.84) 

10(.03) 

(.80.88) 


too  high 
cover 
too  low 
int.val 

161.05) 

2571.86) 

271.09) 

(.82..90) 

341.11) 

2181.73) 

4S1.16) 

(.68.-78) 

111.04) 

264(.88) 

251-08) 

(.84,.92) 

491.16) 

237(.79) 

14(.05) 

(.74,-84) 

251.08) 

2311.77) 

441.15) 

(.72,.S2) 

too  high 
cover 
too  low 
int.val  i 

71.02) 

283(.94) 

101.03) 

(.92.-97) 

161.06) 

2571.86) 

271.09) 

(.822.90) 

41.01) 

2841.95) 

12(.04) 
(.92  .97) 

22(  .07) 
269  .90) 
9(.03) 
(.86, .93) 

81.03) 
2641.8S) 
2S  .09) 
1.S4..92) 

161.05) 
245(.S2 
39  .13) 
(.77..S6) 

51.02) 

2831.94 

12(.04 

(.92.97 


161.05) 

225(.75) 

59(.20) 

(.70,.80) 

161.05) 

2451.82) 

(.771.86) 


411.13) 
2431.81) 
16(.07) 
(.77,. 85) 

20(.07) 
271  .90) 
91.03) 
(.87,. 94) 


23(.OS) 

2311.77) 

46(.15) 

(.72..S2) 

71.02) 
2631.88) 
30(.10) 
(.83  .92) 


too  high 
cover 
too  low' 
int.val 

too  high 
cover 
too  low 
int.val 


201.07) 

257(,86) 
23  .08) 
(.82. .90) 

31.01) 

2971.99) 

01.00) 
(.98. .99) 


201.07) 

1891.63) 
9 1(.  30) 
(.58, .68) 

274i°9l\ 
23(.08) 
1.88, .95) 


581.19) 
2311.77) 
1 1  .04) 
(-72,-82) 

261.09) 
269  .90) 
51.02) 
(.86,  .93) 


39  .13) 
2421.81) 
19(.06) 
(.76,-85) 

191.06) 

2751.92) 

($;!$) 


181.06) 
2261.75) 
56  .19) 
(-70..80) 

32  .11) 
(.S3. .92) 


too  high 
cover 
too  low 
int.val 

26(.09) 

274(.91) 

01/00) 
(.88. .95) 

19^1. 6o) 

981.33) 
(.60,-71 ) 

651.22) 
226  .75) 
91.03) 
(.7l,.S0) 

361.12) 
247  .82) 
171.06) 
(.78, .87) 

151.05) 

229  .76) 
561.19) 
(•72..81) 

too  high 
cover 
too  low 
int.val 

61.02) 
2941.98) 
oi.  oo) 
(.96, .99) 

01.00) 

2021.67) 

98(.33) 
(.62, .73) 

411.14) 
2561.85) 
31.01) 
(.81, .89) 

16(.05) 

2781.93) 

6(.t)2) 

(.90..96) 

6(.02) 

2601. S7) 
341.11) 

(.83, .91)  | 

TABLE  16 


LENGTH  OF  C.I  FOR  In  PfD>  t}(HYPOEXP  MODEL) 
N  -  20.0  -  2/3.PJ  -  2.0, p2  -  10.0, p3  -  2.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J  K 

B.T 

0.5 

80  % 

0.3368 

(0.088) 

0.2706 

(0.088) 

0.1951 

(0.042) 

0.2539 

(0.149) 

0.2382 

(0.105) 

90  % 

0.4186 

(0.110) 

0.3466 

(0.119) 

0.2505 

(0.054) 

0.3306 

(0.194) 

0.3249 

(0.143) 

1.0 

80  % 

0.6008 

(0.158) 

0.5234 

(0.150) 

0.3892 

(0.084) 

0.3918 

(0.161) 

0.3675 

(0.124) 

90  % 

0.7494 

(0.197) 

0.6858 

(0.208) 

0.4998 

(0.108) 

0.5103 

(0.209) 

0.4872 

(0.174) 

1.5 

80  % 

0.9084 

(0.310) 

0.8374 

(0.465) 

0.5840 

(0.129) 

0.6753 

(0.305) 

0.6142 

(0.213) 

90  % 

1.1331 

(0.387) 

1.1701 

(0.834) 

0.7500 

(0.166) 

0.8794 

(0.397) 

0.8144 

(0.29S) 

2.0 

SO  % 

1.3743 

(0.845) 

1.3323 

(1.078) 

0.7S00 

(0.183) 

0.9871 

(0.470) 

0.8909 

(0.316) 

90  % 

1.6939 

(0.930) 

2.1678 

(2.013) 

1.0016 

(0.235) 

1.2855 

(0.612) 

1.1786 

(0.443) 

3.0 

SO  % 

2.5888 

(1.497) 

2.9300 

(2.330) 

1.1776 

(0.357) 

1.6313 

(0.810) 

1.4596 

(0.534) 

90  % 

3.1180 
<  1.550) 

4.'553 

(2.863) 

1.5122 

(0.462) 

2.1245 
( 1.055) 

1.9376 

(0.737) 

4.0 

80  % 

4.2251 

(2.031) 

3.3399 

(2.307) 

1.5894 

(0.763) 

2.2834 

(1.154) 

2.0356 

(0.749) 

90  % 

4.7925 

(1.883) 

4.5249 

(3.430) 

2.0410 

(0.980) 

->  97^9 

1 1.503) 

2.7125 

(1.038) 

IV.  CONCLUSIONS 


This  thesis  considers  the  problem  of  estimating  the  log  probability  for  a  semi- 
Markov  process  which  does  not  enter  a  particular  state  before  time  t. 

Simulation  is  used  to  study  procedures  for  obtaining  confidence  intervals  for  the 
In  probability  a  semi-  Markov  process  enters  a  fixed  state  after  time  t.  The  data  arise 
from  observing  N  individuals.  Three  estimators  and  associated  confidence  interval 
procedures  are  compared:  The  three  estimators  use  different  amounts  of  information 
about  the  process.  The  maximum  likelihood  estimator  and  its  normal  confidence 
interval  uses  the  most  information.  An  estimator  based  on  the  observed  first  passage 
times  uses  the  least.  An  estimator  based  on  an  exponential  approximation  to  the 
survivor  function  of  the  first  passage  time  uses  an  intermediate  amount  of  information: 
confidence  intervals  for  this  last  estimator  are  obtained  using  jackknife  and  bootstrap 
procedures.  The  simulation  results  indicate  the  following. 

(1)  Larger  numbers  of  individuals  resuit  in  smaller  average  confidence  interval 
lengths. 

(2)  The  binomial  confidence  interval  tends  to  overcover.  This  is  true  in  general 
since  the  target  coverage  probability  is  a  lower  bound  on  the  true  coverage 
probability.  They  also  tend  to  have  larger  average  length. 

(3)  The  maximum  likelihood  confidence  intervals  have  the  correct  coverage  if  the 
model  on  which  they  are  based  is  correct.  If  the  model  is  incorrect  they  can 
either  overcover  or  undercover. 

(4)  If  a  jackknife  interval  uoes  not  cover  the  true  value  the  interval  will  tend  to 
be  too  high'. 

(5)  The  confidence  intervals  using  the  asymptotic  renewal  estimator  have  the 
correct  coverage  for  largish  t. 

(6)  The  bootstrap  confidence  interval  requires  much  more  computation  than  the 
jackknife  confidence  interval.  Since  the  results  for  the  jackknife  and  bootstrap 
procedures  are  sirruliar,  ease  of  computation  appears  to  favor  the  jackknife 
procedure. 
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APPENDIX  A 

TRUE  PROBABILITY  TABLE  I 


TABLE  17 

EXPONENTIAL  MODEL  P(D>T1ME}. 


time 

pL  =  1.0  .p,=  10.0 

6  =  0.5 

pj=  1.0  ,p;=  2.0 

0  =  0.5 

Pj  =  1.0.p,=  10.0 

0  =  0.6667 

1.0000 

1.0000 

1.0000 

0.7866 

0.7968 

0.6203 

0.6503 

1.5 

0.4891 

0.5351 

0.3799 

2.0 

0.3857 

0.4415 

0.2753 

2.5 

0.3042 

0.3646 

0.1995 

3.0 

0.2399 

0.3012 

0.1446 

3.5 

0.1891 

0.2488 

0.1048 

4.0 

0.1492 

0.2055 

0.0760 

4.5 

0.1176 

0.1698 

0.0551 

5.0 

0.0928 

0.1403 

0.0399 

5.5 

0.0731 

0.1159 

0.0289 

6.0 

0.0577 

0.0957 

0.0210 

6.5 

0.0455 

0.0791 

0.0152 

7.0 

0.0359 

0.0653 

0.0110 

7.5 

0.0283 

0.0540 

0.0080 

8.0 

0.0223 

0.0446 

0.0058 

8.5 

0.0176 

0.0368 

0.0042 

9.0 

0.0139 

0.0304 

0.0030 

The  sojourn  time  in  state  1  has  an  exponential  distribution  with  mean  1  p 


sojourn  time  in  state  2  has  an  exponential  distribution  with  mean  1  p,. 


APPENDIX  B 

TRUE  PROBABILITY  TABLE  II 


TABLE  18 

HYPOEXPONENTIAL  MODEL  P(D>TIME}. 


time 

Pl  =  2.0  ,p2=2.0 

pr  =  2.0  ,0  =  0.5 

Pj  =  2.0,  p,  =  10.0 

p3  =  2.0,0  =  0.666' 

0.0 

1.0000 

1.0000 

1.0000 

0.5 

0.8652 

0.8670 

0.8214 

1.0 

0.6743 

0.6897 

0.5788 

1.5 

0.5157 

0.5519 

0.3936 

0.3930 

0.4465 

0.2652 

2.5 

0.2992 

0.3268 

0.1783 

3.0 

0.2278 

0.2952 

0.1197 

3.5 

0.1734 

0.2402 

0.0804 

4.0 

0.1320 

0.1954 

0.0540 

4.5 

0.1005 

0.1590 

0.0363 

5.0 

0.0765 

0.1293 

0.0244 

5.5 

0.0582 

0.1052 

0.0164 

0.0443 

0.0856 

0.0110 

6.5 

0.0337 

0.0697 

0.0074 

7.0 

0.0257 

0.0567 

0.0050 

r.5 

0.0196 

0.0461 

0.0033 

8.0 

0.0149 

0.0375 

0.0022 

8.5 

0.01 13 

0.0305 

0.0015 

9.0 

0.0086 

0.0248 

0.0010 

The  sojourn  time  in  state  1  is  the  sum  of  exponential  distribution  with  mean  I  pt 
and  mean  l.p3.  The  sojourn  time  in  state  2  has  an  exponential  distribution  with  mean 

l:p2. 
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APPENDIX  C 

AVERAGE  LENGTH  OF  C.I  FOR  P{D>T} 


Recorded  in  the  Appendix  are  the  average  length  of  confidence  intervals  for 
P{D>  t}  corresponding  to  those  given  in  the  thesis  for  In  P[D>  t}.  If  a  confidence  limit 
is  greater  than  1,  it  is  set  equal  to  1.  If  a  confidence  limit  is  less  than  0,  it  is  set  equal 
to  0.  The  results  are  similar  to  those  in  the  previous  tables. 


TABLE  18 

LENGTH  OF  C.I  FOR  P!D>t}  (EXPONENTIAL  MODEL) 
N  =  50.  0  =  0.5,  pj  =  1.0.  p2=  10.0 


time 

C.I  (%) 

BIN 

NOR  ' 

MLE 

J.K 

B.T 

0.5 

80  % 

0.1651 

(0.U14) 

0.1469 

(0.014) 

0.0692 

(0.007) 

0.1259 

(0.058) 

0.1044 

(0.027) 

90  % 

0.2054 

(0.018) 

0.1887 

(0.019) 

0.0889 

(0.009) 

0.1669 

(0.077) 

0.1330 

(0.033) 

10 

80  % 

0.1915 

(0.006) 

0.1741 

(0.007) 

0.1078 

(0.007) 

0.1278 

(0.036) 

0.1154 

(0.O15) 

90  % 

0.2389 

(0.008) 

0.2236 

(0.009) 

0.1384 

(0.009) 

0.1694 

(0.048) 

0.1483 

(0.01S) 

1.5 

80  % 

0.1967 

(0.002) 

0.1795 

(0.003) 

0.1264 

(0.005) 

0.1381 

(0.035) 

0.1289 

(0.O14) 

90  % 

0.2455 

(0.003) 

0.2305 

(0.003) 

0.1623 

(0.006) 

0.1830 

(0.047) 

0.1650 

(0.019) 

2.0 

80  % 

0.1913 

(.0.007) 

0.1740 

(0.007) 

0.1320 

(0.002) 

0.1427 

1 0.035) 

0.1334 

(0.O15) 

90  % 

0.2387 

(0.008) 

0.2234 

(0.009) 

0.1695 

(0.002) 

0.1892 

(0.046) 

0.1711 

(0.019) 

3.0 

80  % 

0.1706 

(0.012) 

0.1526 

(0.013) 

0.1220 

(0.008) 

0.1325 

(0.031) 

0.1233 

(0.016) 

90  % 

0.2123 

(0.016) 

0.1960 

(0.016) 

0.1567 

(0.010) 

0.1756 

(0.041) 

0.1569 

(0.021) 

4.0 

80  % 

0.1459 

(0.017) 

0.1271 

(0.017) 

0.1007 

(0.013) 

0.1092 

(0.028) 

0.1009 
<0.0  IS) 

90  °  o 

0.1820 
<0.o21 ) 

0.1631 

<0.023) 

0.1293 

(0.017) 

0.1447 

(0.037) 

0.1283 

(0.023) 

45 


TABLE  20 

LENGTH  OF  C.I  FOR  P(D>  t)  (EXPONENTIAL  MODEL) 
N  -  50,  0-2, 3,  p'  -  1.0,  p2-  10.0 


time 

C.I  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.1719 

(0.012) 

0.1601 

(0.012) 

0.0855 

(0.008) 

0.1873 

(0.074) 

0.1645 

(0.044) 

i 

90  % 

0.2160 

(0.015) 

0.2056 

(0.016) 

0.1098 

(0.010) 

0.2482 

(0.098) 

0.2158 

(0.056) 

1.0 

80  % 

0.1962 

(0.003) 

0.1790 

(0.004) 

0.1223 

(0.006) 

0.1523 

(0.042) 

0.1479 

(0.023) 

90  % 

0.2448 

(0.004) 

0.2298 

(0.005) 

0.1571 

(0.007) 

0.2019 

(0.056) 

0.1897 

(0.031) 

1.5 

SO  % 

0.1911 

(0.007) 

0.1737 

(0.007) 

0.1320 

(0.002) 

0.1542 

(0.043) 

0.1474 

(0.019) 

90  % 

0.2383 

(0.009) 

0.2231 

(0.009) 

0.1694 

(0.002) 

0.2044 

(0.057) 

0.1907 

(0.025) 

:.o 

80  % 

0.1772 

(0.011) 

0.1595 

(0.012) 

0.1269 

(0.006) 

0.1496 

(0.042) 

0.1397 

(0.019) 

90  % 

0.2208 

(0.015) 

0.2048 

(0.016) 

0.1629 

(0.008) 

0.1982 

(0.055) 

0.1803  1 

(0.024) 

' 

3.0 

80  % 

0.1435 

(0.019) 

0.1245 

(0.020) 

0.0996 

(0.013) 

0.1178 

(0.037) 

0.1073 

(0.022) 

90  % 

0.1780 

(0.024) 

0.1598 

(0.026) 

0.1279 

(0.017) 

0.1561 

(0.049) 

0.1372 

(0.028) 

4.0 

80  % 

0.1140 

(0.021) 

0.0932 

(0.024) 

0.0701 

(0.015) 

0.0812 

(0.031) 

0.0733 

(0.021) 

-  .. 

90  % 

0.1409 

(0.026) 

0.1185 

(0.032) 

0.0900 

(0.020) 

0.1076 

(0.041) 

0.0934 

(0.028) 
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TABLE  22 

LENGTH  OF  C.I  FOR  P{D>t)  (HYPOEXP  MODEL) 
N  -  50,6  =  0.5,  pj*  2.0,  pUlO.O,  p3-  2.0 


time 

C.I  {%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  % 

0.1417 

(0.01S) 

0.1217 

(0.019) 

0.0692 

(0.006) 

0.0813 
(0.035  > 

0.0722 

(0.018) 

90  0/o 

0.1758 

(0.022) 

0.1574 

(0.024) 

0.0888 
( 0.008) 

0.1077 

(0.047) 

0.093 1 

(0.022)  ! 

1.0 

SO  % 

0.1855 

(0.009) 

0.1680 

(0.010) 

0.1078 

(0.007) 

0.1  no 

(0.032) 

0.1024 

(0.014) 

90  0/o 

0.2313 

(0.012) 

0.2157 

(0.013) 

0.1384 

(0.009) 

0.1472 

(0.042) 

0.1317 

(0.019) 

1.5 

SO  % 

0.1964 

(0.003) 

0.1792 

(0.003) 

0.1265 

(0.004) 

0.1344 

(0.035) 

0.1245 

(O.Olol 

90  % 

0.2451 

(0.004) 

0.2301 

(0.004) 

0.1625 

(0.006) 

0.1182 

(0.047) 

0.1597 

(0.021) 

2.0 

SO  % 

0.1920 

(0.006) 

0.1747 

(0.006) 

0.1322 

(0.002) 

0.1413 

(0.036) 

0.1308 

(0.017) 

90  % 

0.2395 

(0.008) 

0.2243 

(0.008) 

0.1698 

(0.002) 

0.1872 

(0.048) 

0.1676 

(0.021) 

3.0 

SO  % 

0.1675 

(0.014) 

0.1495 

(0.015) 

0.1223 

(0.007) 

0.1265 

(0.032) 

0.1174 

(0.019) 

90  % 

0.2185 

(0.01S) 

0.1919 

(0.015) 

0.1573 

(0.009) 

0.1677 

(0.043) 

0.1496 

(0.022) 

4.0 

80  % 

0.1390 

(0.018) 

0.1198 

(0.019) 

0.1009 

(0.011) 

0.0988 

(0.028) 

0.0919 

(0.019) 

90  % 

0.1723 

(0.023) 

0.1537 

(0.025) 

0.1296 

(0.015) 

0.1310 

(0.037) 

0.1168 

(0.023)  i 
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TABLE  23 


LENGTH  OF  C.l  FOR  P{D>  t}  (HYPOEXP  MODEL) 
N  -  20.  0  -  0.5,  pj  =  2.0,  p2=  10.0,  p3  -  2.0 


time 

C.l  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

0.5 

80  °  o 

0.2328 

10.044) 

0.1839 

(0.057) 

0.1108 

(0.016) 

0.1260 

(0.051) 

0.1177 

(0.044) 

90  % 

0.2836 

(0.053) 

0.2314 

(0.077) 

0.1423 

(0.021) 

0.1641 

(0.067) 

0.1592 

(0.065) 

1.0 

80  % 

0.3005 

(0.023) 

0.2614 

(0.026) 

0.1710 

(0.016) 

0.1736 

(0.050) 

0.1633 

(0.032) 

90  % 

0.3671 

(0.029) 

0.3354 

(0.034) 

0.2196 

(0.021) 

0.2261 

(0.066) 

0.2105 

(0.043) 

1.5 

80  % 

0.3163 
(0.O1 1) 

0.2785 

(0.011) 

0.1990 

(0.010) 

0.2110 

(0.053) 

0.1954 

(0.031) 

90  % 

0.3869 

(0.013) 

0.3576 

(0.015) 

0.2556 

(0.012) 

0.2747 

(0.069) 

0.2502 

(0.040) 

2.0 

SO  % 

0.3104 

(0.016) 

0.2722 

(0.017) 

0.2067 

(0.006) 

0.2195 
i  0.050) 

0.2017 

(0.032) 

90  % 

0.3796 

(0.020) 

0.3495 

(0.023) 

0.2654 

(0.007) 

0.2859 

(0.065) 

0.2571 

(0.038) 

3.0 

80  % 

0.2709 

(0.035) 

0.2289 

(0.040) 

0.1896 

(0.020) 

0.1917 

.0.047) 

0.1756 

(0.037) 

90  % 

0.3304 

(0.043) 

0.2924 

(0.054) 

0.2434 

(0.025) 

0.2497 

(0.062) 

0.2211 

(0.045) 

4.0 

80  % 

0.2276 

(0.046) 

0.1771 

(0.061) 

0.1560 

(0.029) 

0. 1470 
(0.047) 

0.1354 

(0.040) 

90  % 

0.2774 

(0.055) 

0.2223 

(0.082) 

0.2001 

(0.038) 

0.1914 

(0.061) 

0.1696 

(0.050) 

TABLE  24 


LENGTH  OF  C.l  FOR  P{D>t}  (HYPOEXP  MODEL) 
N«  50,  0-2/3,  p1  =  2.0,  p2=  10.0,  p3  =  2.0 


time 

C.l  (%) 

BIN 

NOR 

MLE 

J.K 

B.T 

80  % 

0.1375 

(0.016) 

0.1349 

(0.050) 

0.1237 

(0.034) 

90  % 

0.1938 

(0.019) 

0.1766 

(0.020) 

0.1788 

(0.066) 

0.1624 

(0.045) 

1.0 

0.1330 

(0.036) 

0.1257 

(0.017) 

0.2274 

(0.006) 

0.1763 

(0.048) 

0.1615 

(0.021) 

1.5 

80  % 

0.1324 

(0.001) 

0.1435 

(0.038) 

0.1346 

(0.014) 

m 

0.2397 

(0.008) 

0.1722 

(0.019) 

2.0 

80  % 

0.1747 

(0.013) 

0.1569 

(0.014) 

0.1351 

(0.035) 

0.1257 

(0.015) 

0.2176 

(0.017) 

0.1791 

(0.046) 

0.1606 

(0.020) 

3.0 

0.1327 

(0.019) 

0.1132 

(0.021) 

0.0998 

(0.011) 

0.0957 

(0.028) 

0.0885 

(0.018) 

90  % 

0.1644 

(0.024) 

0.1451 

(0.028) 

00  281 
(0)14) 

0.1268 

(0.037) 

0.1123 

(0.022) 

4.0 

0.0544 

(0.016) 

90  % 

0.0690 

(0.021) 

51 


TABLE  25 


LENGTH  OF  C.I  FOR  P(D>t}  (HYPOEXP  MODEL/ 
N  -  20,  6  =  2,  3,  2.0,  P2=  10.0,  P3  =  2.0 


ro 
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